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We consider the thermodynamic properties of the quasi-two-dimensional spin-half Heisenberg 
ferromagnet on the stacked square and the stacked kagome lattices by using the spin-rotation- 
invariant Green's function method. We calculate the critical temperature Tc, the uniform static 
susceptibility x, the correlation lengths and the magnetization Al and investigate the short-range 
order above Tc- We find that Tc and M at T > are smaller for the stacked kagome lattice which 
we attribute to frustration effects becoming relevant at finite temperatures. 
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Introduction: Quasi-two-dimensional magnets have at- 
tracted much attention in recent yearsi^ According to 
the Mermin- Wagner theorem,'^ strictly two-dimensional 
(2D) Heisenberg magnets do not possess magnetic long- 
range order (LRO) at any finite temperature. Therefore, 
the exchange coupling between planes is crucial for the 
existence of a finite critical temperature Tc- Though 
most of the quasi-2D magnetic insulators are antiferro- 
magnets there are also some quasi-2D ferromagnetic insu- 
lators hke K2CUF4, La^BaCuO^, Cs2AgF4m> The calcu- 
lation of thermodynamic properties for quasi-2D magnets 
is an important issue (S^SiLSiS in particular with respect to 
the interpretation of experimental results. Especially in 
the limit of weak interlayer coupling the evaluation of 
Tc is challenging. Recently accurate values for the crit- 
ical temperature of a Heisenberg magnet on a stacked 
square lattice by Monte-Carlo calculationgiS and on the 
cubic lattice by high-temperature series expansiorJi have 
been presented. Another promising approach to cal- 
culate the thermodynamics of quasi-2D magnetic sys- 
tems is a spin-rotation invariant Green's function method 
(RGM)i^'^d3,i4,i5,i6,i7,i8,i9,20 rj.^^ j^q^ ^^i^^^ Consis- 
tent description of LRO as well as of short-range order 
(SRO) at arbitrary temperatures in magnetic systems of 
arbitrary dimension. 

In this paper we use the ROM for the calculation 
of thermodynamic quantities of the quasi-2D spin-half 
Heisenberg ferromagnet on the stacked square and the 
stacked kagome lattice. Both lattices have the same co- 
ordination number, but differ in their lattice geometry 
within the layers. Note that corresponding RGM treat- 
ments for the respective stacked antiferromagnets were 
given in Refs. ilQ and The Heisenberg antiferro- 

magnet on the kagome lattice is a canonical model to 
study frustration effects in 2D spin-half systems, see, e.g. 
Refs. XJ:''' 21 22 23. Due to strong frustration there is 
most likely no magnetic LRO in the kagome antiferro- 
magnet. It has been argued recentlyi^ that even the 




FIG. 1: Sketch of one kagome layer in the congruently 
stacked lattice with the in-plane geometric basis vectors 
ai = (0,2,0), 32 = (\/3, 1,0). The out-of-plane basis vector 
as = (0, 0, 1) is not shown. The geometrical unit cell contains 
three sites (spins) labeled by a running index a = 1,2, 3. 

increase of the dimension in a stacked kagome antifer- 
romagnet does not lead to magnetic LRO. 

In this paper we consider layered ferromagnets. Then 
geometric frustration is irrelevant in the ground state. 
However, at finite temperatures excited states with an- 
tiferromagnetic correlations contribute to the partition 
function. The energy of such states is influenced by 
triangular spin configurations and thus the frustrating 
geometry of the kagome lattice may become more and 
more relevant with increasing temperature. By compar- 
ison with the unfrustrated stacked square lattice we can 
discuss these frustration effects. 

We consider the Heisenberg model of N spins 1/2 

where the sum runs over all J\f unit cells (labeled by m 
and n) and all spins within a unit cell (labeled by running 
indices a and (3, see Fig. ^for the kagome case). The ex- 
change coupling Jma,nf3 is uonzcro for nearest neighbors, 
only. The ferromagnetic exchange parameter within the 
layers J|| is fixed to J|| = — 1. 

Spin- Rotation- Invariant Green's Function Method: 
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The rotation-invariant decouphng scheme introduced 
by Kondo and YamajiiS goes one step beyond the 
random-phase approximation (RPA). An advantage of 
this method is the possibiUty to treat magnetic order- 
disorder transitions driven by quantum fluctuations as 
wefl as by thermal fluctuations and to provide a good de- 
scription of magnetic SRO,^^i^'^i^4i^''^i^"i'7ii>^ii V" In this 
paper we only illustrate some basic features of the RGM. 
For more details the reader is referred to Refs. mil 
and in particlar to Refs. [T^IT^l20l where the 
corresponding stacked antiferomagnets are considered. 

To evaluate all relevant correlation functions {Sma^n/s) 
and several thermodynamic quantities we have to cal- 
culate a set of Fourier-transformed commutator Green's 
bmctions {(S^^; S^p))^^ which are connected with the dy- 
namic spin susceptibilities by Xq~^(w) = -(('S'+„; S~ij))uj- 
Using the equations of motion and supposing spin- 
rotational invariance, i.e. {S^a) — 0' obtain 

^'^ii^qa'^ ^qp))'-' — ii'^^qa^ ^qp]) + ( ( ~ '^'qa ! 'S'q/3 ) ) w • 

The operator —S^^ = [[S^^, H], H] written in site 
representation contains products of three spin operators 
along nearest-neighbor sequences which are treated 
in the spirit of the decoupling scheme of Shimahara 
and Takadapi^ For example, the product S^S'l^S^ is 
replaced by riA,B(SJiSg)S^ + r]A,c{S^Sc)SB, where 
A, B, C represent spin sites. The vertex parameters 77^^^ 
are introduced to improve the approximation. In the 
minimal version of the RGM we introduce just as many 
vertex parameters as independent conditions for them 
can be formulated, i.e. we have and ri± attached 
to intralayer and interlayer correlators, respectively. 
After the decoupling we can write the equation of 
motion in a compact matrix form omitting the running 
indices a and /3, (w^ — Fq) Xq~ {^) = —Mq, where Fq 
and Mq are the frequency and momentum matrices, 
respectively. Their dimension equals the number of 
spins in the respective unit cell. The solution for 
Xq~' (lo) can be formulated in terms of the common 
set of normalized eigenvectors \jq} of Fq and Mq and 
reads Xq"M = "ELi ^T^^i |jq) (jq|- The case 
J|| > for the stacked kagome lattice has been discussed 
previously^^ in detail. Hence, it is not necessary to 
present the lengthy expressions for the eigenvalues m^q 
of Mq and w|q of Fq once again. The corresponding 
expressions for rrijq and a;|q in case of the stacked 
square lattice are even simpler, since only one eigen- 
vector |lq) exists. For both lattices the eigenvalues 
rrijq and uj'jq contain eight parameters which must 
be determined selfconsistently, namely the spin-spin 
correlators ci,o,o, ci,i,o, C2,o,o, co,o,i' ci,o,i; co,o,2 and the 
two vertex parameters ryy, and r]±. The correlators 

are defined by Ck,i,m = cr = (•S'qS'j^) = 2(SoSr)/3, 
where R — khi + lh2 + mba (the are the basis 
vectors of the cubic lattice) for the stacked square 
lattice, but R = kai/2 + /a2/2 -|- ma^ for the stacked 
kagome lattice, cf. Fig. ^ The spectral theorem24 
for the correlators Ck,i,m yields seven equations in- 



cludin g the sum rule Co,o,o = 1/2. According to 
Refs. II6II8I and an eighth equation is obtained 
requiring that the matrix of the static susceptibility 
Xq~ — Xq~ ('^ — 0) has to be isotropic in the limit 
q ^ 0, i.e., lim,^(^,_o X^"|,^^o = 1™?^"^° ^^L^^^^^o' 
In the RGM scheme magnetic LRO is detected 
,3^13,14,15,16,19 iim|R|^^ (SoSr) ^ 0, i.e. the correlation 
function contains a long-range part Cqq/j, introduced by 
Shimahara and Takada and called condensation parti^, 
leading to (SmaSn^) = 7^ Eq^q •S'qQ/ge"*'^'"""'"'' + 
|C'Qa/3e~'*^'"'"°'"'', where Q is the magnetic wave vector 

and Sq^p = I IS^ (1 + T^j:^) (« Ijq) O'q 

The existence of nonzero Cqq,^ is accompanied by a 
diverging static susceptibility Xq~ at q = Q. This 
gives an additional equation for one Cqq^. If the 
unit cell contains more than one spin, the further 
relations for the various condensation terms corre- 
sponding to different Green's functions are given bj*2Si 
Cqa/s/Cq^a = limq-qXqa/s/Xq^^- condensation 
term CqoP is related to the order parameter M^^ by 
M^^ = 3|Cqa/3|/2 with M^^ being the square of the 
magnetization^^ of a certain site a in the unit cell. In 
a short-range ordered phase the condensation terms are 
zero. For temperatures lower than the respective Curie 
temperature Tc both systems exhibit ferromagnetic LRO 
with Q — (0,0,0). In the kagome case we have three 
spins per unit cell, and in principle six nonequivalent 
Green's functions have to be considered. However, due 
to symmetry we have Cqa^/Cq-y^ = 1, i.e. all three sites 
in the unit cell carry the same magnetization M = Maa, 
a = 1,2,3. The condensation term is obtained by de- 
manding that the static susceptibility in the long-range 
ordered phase diverges at the magnetic wave vector 
Q, which yields J|| (l - 2ri\\ (4ci_o,o - ci,i,o - C2,o,o)) + 
4J_L (77_lci,o,i — '7||Ci_o,o) = for the stacked kagome 
lattice and J|| (l - 27?[| (5ci^o,o - 2ci4,o - C2,o,o)) + 
4J_L (77_lci,o,i — '7||Ci_o,o) = for the stacked square 
lattice. At zero temperature the solution of our 
set of equations becomes exact, namely wc find 
Ck,i.m = C'qa/s = 1/6 and ryu = rij_ ^ 3/2, and 
consequently (SqSr) = 1/4 and Maa = 1/2. 

Further quantities of interest are the uniform static 
susceptibility x — liniq^oXq /2 and the intralayer and 
interlayer correlation lengths and which can be 
calculated as described in Refs. 13 15 16 18 20. 

Results: The critical temperature in dependence on 
the interlayer coupling J± is plotted in Fig. [21 Note 
that the RGM fulfills the Mermin- Wagner theorem.'^, i.e. 
Tc{Ji_ = 0) = 0. In the inset of Fig. H we show M(T) 
curves for J± = —0.2. Obviously, the magnetization M 
at < T < Tc and the value of Tc for the stacked 
kagome lattice are always smaller than M and Tc for the 
stacked square lattice. The relation ^'^"■some ^ rpsquare 
holds for any value of J± (see Fig. 12). To give a possi- 
ble interpretation for that we have to look at the energy 
spectrum for the strictly 2D square and kagome lattices. 
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FIG. 2: Curie temperature Tc as a function of the interlayer 
coupling J±/J\\ ■ For comparison we give the empirical curves 
according to Eq. (0. The inset shows the magnetization M 
in dependence on temperature T for J±/J\\ = 0.2. 



The ground-state energy (lower bound of the spectrum) 
is Emin = — 0.5A^ for both lattices. However, the upper 
bound for the kagome lattice E,nax — 0.434A^ is much 
lower than the corresponding value E,nax — 0.670A^ for 
the square lattice (see, e.g., Ref. Ujj). This upper bound 
is equal to the ground-state energy of the kagome an- 
tiferromagnet (i.e. J|| = -1-1) multiplied by minus one, 
and its small value signals strong frustration effects being 
relevant for eigenstates with antiferromagnetic spin-spin 
correlations. Hence, one can expect that excited states 
with antiferromagnetic spin correlations have lower en- 
ergy for the kagome ferromagnet resulting in a larger 
contribution to the partition function at T > in com- 
parison with the square-lattice ferromagnet. The obser- 
vation y^'^s"™'^ < j,sqiiare QQj-j-ggpQjjfjg data for 

the susceptibility x ^-nd the correlation lengths i.e., 
X and for the stacked kagome lattice are smaller than 
the corresponding values for the stacked square lattice 
(see below. Fig. 

In the limit J^_/J\\ <C 1 we derive more explicit 
expressions for Tc following Ref. 0. We find Tc — 
— J||7r/(2 ln(— 20rc/9Jj_)) for the stacked square lattice 
and Tc = -3J||7r/(81n(-20Tc/9J_L)) for the stacked 
kagome lattice. Note that these expressions are of similar 
form as corresponding formulas for the Neel temperature 
Tm obtained by selfconsistent spin- wave theoriesASiS We 
compare our data for Tc in Tab. with results for the 
Neel temperature Tm of the stacked square-lattice an- 
tiferromagnet. For the simple cubic lattice (i.e. J_l = 
J||) our result for Tc exceeds the best available value 
for Tc obtained by high-order high-temperature series 
expansion-'^^ by 9% which can be considered as a test 
of the accuracy of the RGM results. This accuracy is 
much better than that of the Schwinger boson approach 
(compare 4th and 5th columns in Tab.Hl). Note that the 
accurate results of Ref. ^ show that T^ > Tc for unfrus- 
trated spin-half Heisenberg magnets, whereas an RPA 



FIG. 3: Inverse intralayer (^y ^) and interlayer corre- 
lation lengths vs. T for J±/J\\ = 0.2. The inset shows the 
inverse uniform static susceptibility x for J±/J\\ =0.2. 
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TABLE I: Curie temperatures Tc / {\J\\\s (s + 1)) of the s = 
1/2 ferromagnet on the stacked square and stacked kagome 
lattice calculated by the RGM. For comparison we show 
some results for the Neel temperature Tjv/ ( J||s(s -I- 1)) of 
the stacked square-lattice antiferromagnet obtained by other 
approaches. For the simple cubic lattice high-order high- 
temperature series expansion— yields Tc / {\J\\\s {s + 1)) — 
1.119andrjv/(|J|||s(s-H)) = 1.259. 

treatmentSi yields Tn — Tc- Therefore, the relation be- 
tween Tc and Tjy can be used as a further test of the 
RGM. Indeed, using the same set of vertex parameters'^ 
for J± = J|| we obtain Tn = 1.1225rc which is in excel- 
lent agreement with Ref. 11. 

The overall dependence of Tc versus J± resembles the 
results of Yasuda et almS, for T/v, see Tab. Therefore, 
we adopt their empirical formula 



|j|il B-Hjjj^^y 

The best agreement with our data is obtained for A — 
2.414 (2.049) and B = 2.506 (2.315) for the stacked 
square (kagome) ferromagnet, cf. Fig. |2] 

The correlation lengths and the uniform static sus- 
ceptibility are shown as functions of T in Fig. |21 for 
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R=( 1,0,0) 
R=( 1,1,0) 
R=(2,0,0) 
R=(0,0,1) 
R=( 1,0,1) 
R=(0,0,2) 




FIG. 4: Short-range intra- and interlayer correlation functions 
(SqSr) for the stacked square lattice with J±/ J\\ = 0.05. 



J^_|J\\ = 0.2. As expected, those quantities diverge at 
Tc and they go to zero for T — > oo. The inverse sus- 
ceptibihty depends linearly on T for temperatures well 
above Tc, but it deviates from the linear Curie- Weifi law 
near Tc (see inset of Fig. 13). 

Finally, we discuss the short-range spin-spin correla- 
tors. As an example we present their temperature de- 
pendence for the stacked square lattice for a weak in- 



terlayer coupling J± = —0.05 in Fig. 0] For this value 
of J_L we have Tc = 0.4388, i.e. for T > 0.4388 the 
long-range part of the spin-spin correlation function van- 
ishes. It is obvious that the short-range intralayer and 
interlayer correlators behave differently. The interlayer 
correlators become very small aXT — Tc whereas the in- 
tralayer spin-spin-correlations are still well-pronounced. 
This behavior corresponds to the data for the intra- and 
interlayer correlation lengths, see Fig. |31 and indicates 
the crossover from three-dimensional magnetic LRO at 
T < Tc to 2D SRO at T > Tc, if |J±| < |J|||. 

Summary: In this paper the thermodynamics of the 
spin-half Heisenberg ferromagnet on the stacked square 
and kagome lattices has been investigated applying a 
spin-rotation-invariant Green's function method. We 
have calculated the Curie temperatures Tc in dependence 
on the interlayer coupling J_l, and simple empirical for- 
mulas for Tc{J±) are obtained from the numerical data. 
Studying short-range spin-spin correlation functions we 
see clearly a dimensional crossover at T ~ Tc for strongly 
anisotropic systems. Comparing the values of Tc, of the 
magnetization, the correlation lengths and of the uni- 
form static susceptibility we come to the conclusion that 
at T > frustration becomes relevant for the stacked 
kagome lattice leading to a weakening of magnetic order 
at finite temperatures and to j^^'^sorne ^ rpsquare^ 
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The stacked square-lattice Heisenberg antiferromagnet was 
treated within the RGM in Ref. However, in that pa- 
per an additional vertex parameter was introduced which 



results in slightly different values for Tjv compared to the 
scheme used in the present paper. 



